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^ ' The dynamics of high partonic density QCD is presented considering, in 

1—5 ! 
C^ I the double logarithm approximation, the parton recombination mechanism 

built in the AGL formalism, developed including unitarity corrections for the 
nucleon as well for nucleus. It is shown that these corrections are under the- 
oretical control. The resulting non linear evolution equation is solved in the 
asymptotic regime, and a comprehensive phenomenology concerning Deep In- 
elastic Scattering like F2, Fl, F^. dF2/dlnQ'^, dF-^/dlnQ"^, etc, is presented. 
The connection of our formalism with the DGLAP and BFKL dynam- 
ics, and with other perturbative (K) and non-perturbative (MV-JKLW) ap- 
proaches is analised in detail. The phenomena of saturation due to shadowing 
;h ' corrections and the relevance of this effect in ion physics and heavy quark 

production is emphasized. The implications to e-RHIC, HERA- A, and LHC 
physics and some open questions are mentioned. 
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I. INTRODUCTION 

The dynamics of the high density Quantum Chromodynamics (hdQCD) is one of the 
present most challeging open questions in high energy physics. The intense theoretical 
and experimental activity towards the understanding of small x (small fraction of proton 
momentum carried by the struck parton) QCD takes place from Deep Inelastic Scattering 
(DIS) at HERA |jT|] to heavy ions collisions (HlC) at RHIC [Q]. This kinematical regime will 
also be tested at LHC in a near future |^. 

Important contribution to the interest of the field is due to the puzzling result obtained by 
HERA at small-x {x < 10^^) @] for the proton structure function F2{x,Q'^). This function 
was observed to increase dramatically as x gets smaller (Fig. 1). In the region of moderate 
Bjorken x {x > 10^^) the Operator Product Expansion (OPE) methods as well as the 
Renormalization Group Equations (RGE) have been applied successfully [^. The evolution 
of quark and gluon distribution functions given by the DGLAP equations is based on 
the summing of the leading powers of a^lnQ^ ~ 1, asln(l/x) << 1, ag << 1, where 
as is the strong coupling constant. The leading ln(l/x) contributions is the case for the 
BFKL equation 0. The procedure known as the double leading logarithmic approximation 
(DLA) corresponds in axial gauges to generate the logarithms by ladder diagrams, whose 
emitted gluons have strongly ordered transverse and longitudinal momenta, summing the 
logs Q;sln(5^1n(l/x). It was shown that the DLA is a common limit between the linear 
dynamics P]. 

The increasing of the parton densities requires a formulation of the QCD at high partonic 
density, where unitarity corrections (UC), not considered in the previous dynamics already 
mentioned, are properly taken into account. In this sense, the small x region, where the 
gluon distribution sets the behavior of the main observables, provides the interface between 
perturbative and non perturbative QCD, or in other words, between hard and soft physics. 
Clearly, both experimentalists and theoreticians are challenged to desentangle, measure 
and formulate the dynamical collective effects that are subjacent to the observed result 



of increasing F2 and the cross section atot at DIS, as x gets smaller [Q. Both evolution 
equations, DGLAP (evolution in InQ^) and BFKL (evolution in ln(l/a;)) as representatives 
of linear dynamics, need control in order to restore unitarity, since the Froissart limit requires 
atot<Cte\n^s §. 

A comprehensive treatment should envolve both linear and non linear regimes. The 
main attempts to develop a formalism for hdQCD are the approaches of Mc.Lerran and 
collaborators (MVJKLW) |jlO|, by Kovchegov (K) |ll| and by Ayala, Gay Ducati and 
Levin (AGL) |]l2|,|l3l. Derived independently, the three methods obtain non linear evolution 
equations for the gluon distribution, at the small x region, describing the onset of hdQCD, 
although considering different degrees of freedom. 

In what follows I will present an introductory review of the subject of hdQCD, the main 
aspects of the formulations to the subject, the connections among them in the asymptotic 
region (x — ;> 0), present the state of art of the phenomenology of small x physics and address 
some open questions. 

II. THE EVOLUTION EQUATIONS 

It will be briefly presented the DGLAP and BFKL dynamics and their predictions for 
the small x regime at DIS. The DIS is the process of interaction of a lepton and a nucleon 
exchanging an electroweak boson producing many particles at the final state, which is a 
hadronic state X. The process is 

l{k) N{p) ^ l'{k') X , (1) 

where k, k', p and p' are the fourmomenta of the initial and final lepton, incident nucleon 
and final hadronic system, respectively. The main variables for this process are Q^ = —q^ = 
— {k — k'Y, which is the square of the transfered momentum, s = {k + p)^, the square of 
the center of mass energy, W"^ = {q + p)"^ = {p'Y'i the square of the center of mass energy 
of the virtual boson-nucleon system. The hard scale is given by q^ (< 0), corresponding to 



the process resolution, and x = Q'^/2p.q, meaning the virtual boson resolves the hadronic 
structure, or the partons, since Aa; ~ l/y/Cp . Are also useful the variables y = q.p/k.p, 
measuring the process inelasticity and u = q.p/mi^, the energy of the virtual boson once 
taken in the target rest frame. 

In terms of the partonic content of the nucleon the structure function, which reflets its 
overall distribution, is given by 

F,{x,Q')=Y.e}[q^ix,Q')], (2) 

/ 

where the sum is over flavours weighted by the respective squared charges (cj). It is this 
function that is object of main experimental studies, specially at HERA, at the small longi- 
tudinal nucleon fraction of momentum, or small x (see Fig. [^) |1l4| . 
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FIG. 1. 



The quark distribution function can be shown to evolve as 
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where P, 



qq 



C, 



' 1-2 



(with self-energy corrections over the k propagator) is one of the 



splitting functions Rij ( Pgg = Cp '^^ ~^' , etc), describing the transition between the quark 
state i to the quark state j, from fraction of momentum Xi to x. The above evolution refers 
to the non singlet quarks distribution where sea quarks and gluon distribution are uncoupled 



qNs{x, Q^) = qi{x, Q^) - qj{x, q^ 



(4) 



The singlet quark distribution is given by 



f 



(5) 



where the gluon distribution is coupled to the quark one. Now the evolution equations, in 
the linear regime, read for the quarks 
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and for the gluon distribution 



dg{x,Q^) as{Q' 
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The Eqs. (^ and (0) were independently derived by Dokshitzer, Gribov and Lipatov, and 
by Altarelli and Parisi, known as DGLAP equations in leading order. The perturbative QCD 
evolution, in the linear sector is governed by DGLAP equations, moreover a suitable non 
perturbative inicial condition, extracted from the experiment for a given boson virtuality. 

It can be shown by sucessive derivations that qNs{x,B) ~ J^ni'^s^)"', £ = InQ^, which 
corresponds to the emission of n gluons, showing that the DGLAP equations resum the lead- 
ing InQ^. This can be understood as ladder diagrams with a strong ordering in transverse 
momenta k±, i.e., Ql « kf « ... « k^ « Q^. The scale Ql is the cut, or transition 



value between perturbative and non perturbative physics. It was shown by Gribov |T5| that 
this result is gauge independent once one considers the leading logarithm approximation. 

At small-x the gluons dominate, since Pg^\z) ~ ^^, and the parton distributions 
have the general behavior xpi{x.,Q'^) ~ x^^, A > 0. More likely for initial condition 
xpi(x, Qo) ~ Const and xpi{x,Q'^) ~ expJ\n{\nQ^)\nl/x, known as double leading log- 



arithm approximation (DLA). From that is clear that DGLAP predicts the increase of the 
gluon distribution function, and of the structure function F2 with the decreasing of x, whose 
relation in this kinematical regime is given by 

being equal to ^xg{x, Q^) for Uf = 3. The DLA implies strong ordering in x and fcr, i- 
e., xi » X2 » .... » Xi-i » Xi » X and kr^ « kx^ « •••^Ti-i << kxi « Q^, the 
resum of logs of a^ ln(l/a;) InQ^, having as region of validity a^ << 1, a^ ln(l/2;) << 1 and 
a,ln(l/x)lnQ2 ~ 1. 

The resum of all leading logarithms of Bjorken x, or the energy, is characteristic of the 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation. For very low x values the Ins becomes 
large and as In 1/x ~ 1 and the DLA is not valid. The BFKL evolution equation is proposed 



for an unintegrated gluon distribution function in the transverse momentum variable. Its 
solution grows as a power of the center of mass energy s with the consequent violation of 
the unitarity bound B at very high energies. The cure for this problem was not reached in 



the next to leading order calculation [jT6[, and is still under research for instance, through 
the resumming of all BFKL Pomeron exchanges |T^, for the cross section as well as for the 
structure function. 

The amplitude for the scattering quark-quark with one gluon exchange in the t channel 
at lowest order is given by Ao{s, t) ~ s/t, and for the next order in as the leading terms are 
given by Ins, resulting Ai{s,t) ~ Ao{s,t) Ins. Once one goes to higher orders the number 
of contributing diagrams increases and the calculation gains enormously in complexity ||18|| , 
and the usual procedure is to introduce an effective vertex (Lipatov vertex). It results 
that the general term is An{s,t) ~ Ao{s,t)e"'(t)\n'^{s)/n, where e(t) is a suitable function 
to take care of infrared divergencies, docile under regularization, for instance, dimensional 
regularization. 

Clearly the BFKL evolution is summing the terms a" ln"'(s), where lower order logarithms 
are neglected. The result for the amphtude is 

A{s, t) = Aois, t) f !!W]!!!(£) ^ ,«W , (9) 

n=0 ^ 

with a{t) = 1 + e(t). In this case, there is still the e(t) infrared divergencies to be cured. 

When just the singlet contribution in the t-channel is considered in the BFKL formalism, 
meaning that only Pomeron exchange diagrams are taken into account, the amplitude is 
given by 

ImA{s,t) Q f j2i j2i ^ n \ 

^B{k2,q), (10) 



27r2 
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kl{k-qf 

where Q is the color factor for the process and q is the transfered momentum in the t-channel; 
the functions $« are the impact factors setting the coupling of F to the external particles 
and finally the function F is the perturbative gluon ladder. At leading order it consists 



into the exchange of two gluons but the sum of all terms implies an integral equation for F, 
that is infrared finite for a reggeized gluon ladder. This behavior of the kernel of the BFKL 
equation is connected with the QCD Pomeron and it is the resum of the leading logarithms 
Ins. 

The solution of the BFKL equation predicts the steep growth of the gluon distribution 
with decreasing x as well as the diffusion of the transverse momenta. As well as DGLAP 
equations, the BFKL equation predicts the growing of the cross section in the small-x regime 
since the dynamics of this observable is related with the gluon distribution function. 

From this very brief discussion on the main issues of the linear formalisms for the dynam- 
ics of the parton distributions, it gets clear the need of formal improving in order to include 
the unitarity corrections preserving the Froissart limit. This important aspect of high energy 



physics was pointed out many years ago by Gribov, Levin and Ryskin (GLR) in Ref. [|T^]. I 
will present in the following the main attempts developed in the recent years towards a non 
linear dynamics for high density QCD, as well as the high energy phenomenology provided. 

III. THE QUESTION 

The main question that is addressed once treating hdQCD is how to analytically sepa- 
rate small and large distance contributions to high energy amplitudes in a properly gauge 
invariant formalism. This corresponds to establish the hard and soft scales for the process 
of interest and develop the physical meaningful method to introduce the unitary corrections 
(UC) into the parton dynamics. Once the theoretical need for UC is established we should 
look for their signatures analysing different observables. Besides comparing the predictions 
of the distinct formalisms it is required a common limit between them, probably to be set by 
a saturation scale, (5|. There exists mainly two non linear perturbative approaches [pn]-p!3[] 
and a non perturbative one |Ty]. Although some progress has been made towards their 



connection there is no common analytic solution for the gluon distribution g{x,Q'^) for all 
kinematical range. 
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The physics of hdQCD shakes the parton model and the cherished concept of incoherence 
that is behind the standard calculations. It seems that in order to control the increasing 
of the gluon function some gluon recombination mechanism has to take place as the energy 
gets higher and higher. 

This is a good point to remind that the analysis of the structure functions has already 
given us some surprises, and the previous important one was the EMC effect [^, that has 
as a main result F2 / AF2 7^ 1. This difference is not predicted if one requires complete inco- 
herence of the partons, and reveals the presence of nuclear effects in the structure functions 
where they were not expected. A large literature is devoted to this phenomena, but it is 
interesting to point out that the shadowing behavior noticed in J /'^ production could be 



nicely described [^, as well as the comparison with Drell-Yan processes p2|, for the first 



data at small-x, considering the recombination approach developed by Mueller and Qiu ||2^ 
which is based on the GLR proposals. 

Two main aspects are in order, the control of the increasing of the gluon distribution 
function as an unitarity imperative and the appearance of nuclear effects in high energy 
processes. If this aspect is relevant for fixed target quarkonia production, it is strongly 
important for physics of HERA-A, RHIC and LHC with nuclei. 

IV. THE HIGH DENSITY QCD APPROACHES 

The leading logarithm approximation DLA (related to DGLAP) and LL(l/x) (related 
to BFKL) result both into linear evolution equations for the gluon distribution function. 
The effect of summing large logs in high energy regime implies the increase of the gluon 
distribution function g{x,Q'^) as well as the cross section once x decreases. However, this 
result violates the unitarity of the scattering matrix, a main theorem of relativistic Quantum 
Field Theory, the Froissart theorem 0, which states the cross section cannot grow faster 
than In s. Translated to the DIS, this implies increasing restrictions to the structure function 
and/or total cross section, say lower than In (1/x) and provides a limitation in the x range 
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to the application of linear evolutions in order to get suitable results. 

Intuitively we can associate xg{x,Q'^) to the number of gluons into the nucleon, n^, 
per rapidity unity, y = ln(l/x), with transverse size of order 1/Q. In the hadron-nucleon 
interaction it is the virtual gluon that probes the nucleon structure, in analogy with the 
eletroweak boson in DIS. The virtual gluon- nucleon cross section is 

aG*N{x,Q^) = aoxg{x,Q^) , (11) 

where ctq = cra^g^x = Cte ^a , is the total cross section of the virtual gluon, with virtuality 
Q^, and nucleon gluon interaction. Assuming a^ = nR'jj^^, then aoxg{x,Q^) corresponds 
to the area occupied by the gluons in a nucleon. As x -^ 0, this transverse area may 
be comparable, or even bigger, than t^R\adi following DGLAP or BFKL predictions for 
small X or small Q^. Approaching this regime the gluons may begin to superpose spatially 
in the transverse direction and to interact, behaving not anymore as free partons. These 
interactions should slower, or even stop, the intense growing of the cross section, fixing the 
limit TfRfjj^^ in the small x regime. 

Introducing the function k, with probabilistic interpretation 

n = ^0—;^^^ , (12) 

it is possible to estimate in which kinematical region one can expect modifications in the 
usual evolution equations. So to say, for k, << 1 the system stays at x and Q^ where the 
usual evolution equations (linear) are applicable, governed by individual partonic cascades, 
without interactions among the cascades. 

As K, ^ as, partons from distinct cascades begin to interact due to spatial superposition. 
This specific kinematical regime or the onset of the recombination mechanism was first stud- 
ied by Gribov, Levin and Ryskin |T^ almost twenty years ago, proposing the introduction 



of non linear terms into the evolution equation. 

The region of k — ^ 1 was addressed more recently [T^JT^ (see Fig. |^) and experienced 



considerable development on the theoretical side 0,|TT|], also motivated by HERA results 
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and the great interest on RHIC and LHC future data. This is the kinematical regime that 
requires the QCD dynamics for high partons density. Although the couphng constant a^ is 
still small, allowing in principle the use of perturbative methods, the system is so dense that 
manifestation of non-linear effects are expected, and they are required to be considered in a 
complete formalism. 

The region of k -^ 1 corresponds to partons in a non-equilibrium state and new methods 
are in order to treat the collective phenomena. 



Contour for k = ote for Nuoleon 

r' = 5 GeV' 



20.0 
ln(1/x) 











K =0.6 






= 0.8 






=1.0 

=1.2 


- 




=1.4 


^^ 




HERA 96 


.-><^;- 


~~ — ~~— ~- <^ ^' ''■■'^Z''''^ 






H ERA ^ -X-^t>i^^^^^ 




^^0^^ 


-^^>^ 



-1.0 0.0 1.0 2.0 3.0 4.0 



FIG. 2. 



A. The GLR Formulation 



Gribov, Levin and Ryskin [|19[ introduced the mechanism of parton recombination in 
perturbative QCD for high density systems, expressing this as unitarity corrections included 
in a new evolution equation known as GLR equation. In terms of diagrams it considers the 
dominant non-ladder contributions, or multi-ladder graphs, also denoted fan diagrams. 

The standard QCD evolution is represented by a cascade of partonic decays in the nu- 
cleon. The photon interacts with a parton with fraction of momentum x and virtuality Q^, 
which is the last one of a chain where the partons get slower and with bigger virtuality. 
The scale Qq sets the initial virtuality and at the same time, the limit for perturbative 
QCD applicability, and Q^ is the higher virtuality of the chain. In the transverse plane 

11 



the partons with low fraction of momentum stay in the lower part of the ladder and have 
large transverse size; those with bigger virtuality are on the upper part of the ladder and 
transversally smaller. 

Following DGLAP, the number of partons of low fraction of momentum increases very 
rapidly, which pictorically corresponds to bigger density of individuals in the same allowed 
area, in contrast with a more diluted system at intermediate values of x, far away from the 
possibility of superposition. The transition between these regimes should be characterized 
by a critical value x = xcrit- The same can be argued through BFKL formalism, with the 
difference that in this case the increasing of the partonic distributions, takes place at a fixed 
transverse scale, although the evolution presents the fluctuations in the transverse plane due 
to the characteristic diffusion in BFKL. 

It is important to emphasize that in both linear dynamics only the decay processes 
are considered in the partonic evolution, however we expect that the anihilation mechanism 
should contribute in the low x regime, providing some control of the increasing of the partons 
distribution function. In the linear approach we consider one incident and two emergent 
partons to construct the splitting functions for the decay processes. Now it is the case 
to consider two incident partons and one emergent one, and to express the recombination 
mechanism it is needed a formulation in terms of the probability to recombine two incident 
partons. 

As a first approximation one considers the anihilation probability as proportional to the 
square of the probability to find one incident parton, introducing a non-linear behavior. 

Taking p = ^^^^ as the gluon density in the transverse plane, one has the general 
behavior: for splitting 1—^2, the probability is proportional to tts p, for anihilation 2—^1, 
the probability is proportional to a^p'^/Q^; where l/Q"^ stands for the size of the produced 
parton. For x ^ 0, p increases and the anihilation process becomes relevant. Considering a 
cell of volume Aln(5^Aln(l/x) in the phase space allows one to write the modification of 
the partonic density as 
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where the couphng in the process is given by 7. Expressing in terms of the gluon distribution 
the above equation is 



dlnQ'^dln 1/x 
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(14) 



This equation is the GLR equation |19]. The aheady mentioned work of Mueller and 



Qiu gives 7 = 81/16 for N^ = 3. 

The non-linear corrections correspond to a class of QCD FejTiman diagrams, called fan 
diagrams, formed by a gluon ladder with subsequent subdivisions in gluon ladders, where 
the three ladders vertex is associated with the decay and consists of a sum of several non 
planar diagrams. The overall result carries a minus sign, which is important in order to 
control the growing of the parton distribution once the fan diagrams become relevant, i.e., 
at low X. The lowest part of the diagrams couples to the nucleon and the Eq. (14) resums 
all class of the diagrams represented in Fig. 13. 
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FIG. 3. 

As is clear from Eq. (14), the non-linear term reduces the growing of xg{x, Q'^) at low x, 
in comparison with the linear equations. It is also predicted for the asymptotic region x ^ 
the saturation of the gluon distribution, with a critical line between the perturbative region 
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and saturation region, setting its region of validity (meaning independence of the gluon 
function with the energy). The subject of saturation is very appeahng and there are several 
attempts in the literature today with distinct phenomenological approaches addressing this 



question 24 



In the asymptotic limit one obtains xg{x,Q'^) g^^ = 27^ Q^R^- Since the GLR only 
includes the first non-linear term, although it predicts saturation in the asymptotic regime 
its region of validity does not extend to very high density where higher order terms should 
contribute significantly. 

B. The AGL Formulation 

This approach developed by Ayala, Gay Ducati and Levin (AGL) |T2| , p!3| , intents to ex- 
tend the perturbative treatment of QCD up to the onset of high density partons regime, 
through the calculation of the gluon distribution which is the solution of a non-linear equa- 
tion that resums the multiple exchange of gluon ladders, in double leading logarithm ap- 
proximation (DLA). 

It is based on the development of the Glauber formalism for perturbative QCD ||25|| , 
considering the interaction of the fastest partons of the ladders with the target, nucleon or 
nucleus, since one of the main goals is to obtain the nuclear gluon distribution xg^{x, Q"^)- 
We considered a virtual probe G* that interacts with the target in the rest frame, through 
multiple rescatterings with the nucleons. In this reference frame the virtual probe can be 
interpreted following the decomposition of the Fock states, and its interaction with the target 
occours by the decay of the component gg, as represented in Fig. (|^). 

For small-x this pair has a lifetime much bigger than the nucleus (nucleon) radius and 
the pair is separated by the fixed transverse separation r^ during the interaction, which is 
represented by the exchange of a ladder of gluons strongly ordered in transverse momentum. 

The cross section for this process is given by 



^G*A 



f'dz f^\^r{Q\n,x,z)\'a^^^\ (15) 

JO J IT 
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where G* is a colorless virtual probe with virtuality Q^, z is the probe fraction of energy 
carried by the gluon and '^f* is the wave function of the transversely polarized gluon in the 
probe, a^^^'^^z^r'l) is the cross section of the pair with the target, which was proven for 



perturbative QCD by Mueller in Ref. [25,26 



The lower limit estimation of UC is obtained through the incoherent rescatterings of the 
gluon pair, with the constraint that only the fastest partons of the ladders interact with the 
target. Introducing the transverse impact parameter ht and a profile function for the nucleus 
S(ht) we get 

JQ Jo IT J TT 

2[i-e5<^-'.4A?)^(^0] ^ (16) 

where x' = x/{zrtQ'^), S{bt) may be taken as ^;^e~^*/^^ for a gaussian profile, af^ = ^cr]^, 
where a]^ = ^{^as{4:/rf))7r'^rfxg{x, 4:/rf), and A/rf is a cut for the nonperturbative region. 
For the virtual probe with virtuality Q^ the relation a'^*^{x,Q^) = ( ^q"0 ^5'a(3^) Q^) is 
valid. 

In this approach the gluon pair emission is described in DLA of perturbative QCD and 
from the Feynman diagrams of order a", it should be extracted only the terms that contribute 
with a factor of order (a^ In 1/x In Q'^/Qq)"'- The interaction of the gluon pair with the target 
operates through the exchange of a ladder which satisfies the DGLAP evolution equation in 
the DLA limit. 

It is a working hypothesis that in high energy the successive rescatterings can be taken 
as independent allowing the employ of Glauber formalism, in such a way using the eikonal 
procedure for a relativistic particle propagating in the target. 

Our master equation for the interaction of the gg pair with the target is known as the 
Glauber-Mueller formula and is 

4 /■! dx' r°° (Prt f d\ 



xgA{x, Q ) = — ^ , „ — I / 
n^ Jx x' Ja Q^ Tirf J 



A/Q'^ Txr]: J n 
2[l-a%'{x',4/r^)Sm. (17) 



15 



The pictorial i 



ig- 



(I)- 




FIG. 4. 



Once we perform the impact parameter integration using a gaussian profile function we 
obtain 



xgA{x,Q'^ 



2R\ /•! dx' f^/Ql (fn 



TT^ Jx X' Jl/Q^ TTTJ 



C^\x,{KG{x',rl))^E^{KG{x\ri)) 



(18) 



where C is the Euler constant, E\ is the exponential function and where the kq function 
was introduced as 

3q;.o 



KG{x,r\) 



2R' 



-nrtxg{x,l/rt) 



(19) 



The expansion of Eq. (18) in terms of kq gives as the Born term the DGLAP equation 
in the small x region, the higher order terms corresponding to the unitarity corrections 
naturally implemented in this formalism. 

The estimation of the shadowing effect due to gluon recombination can be immediately 
obtained studing the ratio -Ri = xqa^x, Q^)/Axg'^^^{x, Q"^) presented in Fig. (H), where we 
calculate for Ca and Au, and analysed the behavior of this function in terms of ln(l/x), 
A^/^ and InQ^. We used the GRV parametrization |2^ and adapted the calculation in order 
to have a larger domain of validity in x using 



xg{x, Q^) = xgA[Eq.{18)] + Axg^'^^ix, Q 
-A / / --—xg{x,Q) 

TT Jx JQl Xi Q ■^ 
16 



(20) 



where Axg'^^^^XjQl) is the initial condition. The same procedure could be applied for 
another global parametrization based on DGLAP. 

As expected the UC increase with A, and get smaller as Q^ increases and it is evident the 
importance of the effect as x goes to small values. This allows us to say that the UC should 
be included in the calculations related with the nuclear gluon distribution function, and the 
obtained function xgA{x, Q"^) may be used to set the initial conditions for future experiments. 
For instance in HERA-A, in processes e^A [^ the function xqa could be obtained indirectly 



and employed as an initial condition for the hadronic high energy processes at RHIC and 
LHC. 
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The quarks and gluons distribution were also analysed for the nucleon in this formulation 



|13| , as well as the structure function F2 [^. The motivation for this generalization is the 
availability of HERA data. 

The free interpretation of the Froissart theorem for hadronic processes requires a limit 
for the increasing of the cross section o"^*Ar and F2 with the energy so unitarity is not 
violated. Concentrating the discussion on the k value, k = xg{x,Q'^)(j^^/{Q'^TTR'^) = 



SnasXglx, Q"^) /2Q'^R'^ , which is the probabihty of gluons interactions inside the partonic 
cascade, and R is the radius of the nucleon area occupied by the gluons, we were able to 
obtain R^ = 5 GeV~^, and that k, reaches 1 at HERA, meaning the effects of shadowing 



should be considered in the analysis [|T^,|1^,^]. In the nucleon case, following the same steps 
as before we obtain 

4 /■! dx' f^ (Prt f^ d\ 



TT^ Jx x' J4/Q2 vrr? JO 



4/Q2 vrr^ Jo IT 
2 [1 _ e-i-^'(-''4/r?)s(60] ^ (21) 

and requiring the recovering of DGLAP at DLA we have 

a%'{x, AM) = ^r^xg{x, A/r^) . (22) 

For the quarks, considering the scattering of a virtual photon that decays into a quark- 
antiquark pair, which interacts with the nucleon through the exchange of a ladder we get 

a{j*)=jy^JdW^{z,n)\'aZ^'', (23) 

where "^ is the wavefunction of the qq in the virtual photon [^] . We obtain 

(Ttot = / d\ [1 - e-^^^«(^''^''''*)] . (24) 

Here Qqq is the opacity function that in the Glauber (or eikonal) approach is equivalent 
to il = "" ^°2 ' xg{x,Q'^)S{bt). In doing so we are able to reproduce DGLAP evolution for 
fi < 1, and guarantee the validity of the formulation also for the kinematical region where 

n>i. 

Taking f2 ^ oo and factorizing the bt dependence we obtain the unitarity limit for the 
structure function, having for the InQ^ derivative of F2 , dF2/dlnQ^ < Q"^ R'^/Sn'^. Using 
GRV, the unitarity limit for HERA is reached for Q^ = g2 ^ 1 _ 2 GeV^ {y = In 1/x ~ 9). 
Similarly, for gluons it is Q^ = 1 _ 2 GeV^ {y = In 1/x ~ 7) for HERA [0. 

With the aim to obtain a non-linear evolution equation containing the unitarity correc- 
tions through the inclusion of all the interactions besides the fastest parton from the ladder, 
we differentiate our master equation for the gluon in y = In 1/x and e = InQ^, obtaining 
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dyde =^^1^ + H'^g) + E,{kg)], (25) 

where k^^^^^^x, Q^) = ^§xg'^^^^P{x, Q^) for calculations. 
In terms of k,g the main evolution equation is 

d'^K(y,e) dK,(y,e) NcOg r^ , / x ^ , ^^ ,^„^ 

-^TIT^ + n = ^ P + 1^ '^G + E,{kg)] . 26 

It should be mentioned that large distance effects are absorved in the initial condition for the 
evolution, and situating in a conveninet region of Q^ only short distance effects are present, 
meaning a perturbative calculation is reliable. 



Equations (25) and (26) were derived in Ref. [|12],|13|, and refered for simplicity as AGL 



equation. The main properties of this formulation are: 

• all contributions from the diagrams of order (a^ye)^ are resummed; 

• in the limit k — > the DGLAP evolution in DLA is fully recovered; 

• for K, < 1, and not large, the GLR equation is recovered; 

• for Osye ~ 1 the equation is equivalent to the Glauber formalism. 

The UC are described for the different kinematical regions of k from strictly perturbative 
QCD up to the onset of hdQCD. Non-perturbative effects are not explicitly described and 



this is the object of a distinct formalism MV — JLKM |10] that we will briefly comment in 
a next subsection. In Fig. (^ the comparison between the solutions of the equations AGL, 
GLR, DGLAP and Glauber-Mueller (MOD MF) formula is presented, where the control of 
the growing of the gluon distribution once UC are considered is very evident. 

It was also obtained the asymptotic solution k > 1 of the AGL equation for fixed as 



1^ , P^ as well as for running a^ I^Tj. For high partonic density and y » y^ we obtain 



«:a"(y) = ^ylny^^y. (27) 
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This solution is a good approximation for very small values of x ((9(10 ®)), related with 



THERA physics ||2^, region of a very dense parton system. In terms of the gluon function 
the asymptotic behavior is 



xg{x,Q ) 



2, 2N^Q^R'^ 



3 7r2 



ln(l/a;) 



(28) 



presenting a behavior softer than predicted by DGLAP, meaning a partial saturation. 
For the running a^ the result is [^ 

e 2N^Q^R^ 



xg{x,Q' 



ln(l/x) 



(29) 



1 + e 37r2 

where e = InQ^/Ag^^,. The partial saturation is not modified, and the main difference from 
the previous result occours for small values of e. This confirms the expectation that the UC 



are already relevant before the corrections to leading order pO|,^ . 



C. The Kovchegov Formulation 

The unitarization problem in QCD was addressed as an extension of the dipoles formalism 
for the BFKL equation by Kovchegov [|r^ . This work proposes a non linear generalization of 
BFKL equation, also addressed previously in Ref. P3 by the use of OPE to QCD obtaining 
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the evolution of Wilson line operators. The scattering of a dipole (onium - qq) with the 
nucleon is described by a cascade evolution corresponding to the successive subdivision of 
dipoles from the father dipole. Each dipole has multiple scatterings with the nucleons of 
the target, implying multiple ladders exchange to be resummed in order to obtain the cross 
section of the interaction of the dipole with the nucleus. As a result it is derived the evolution 
equation having the unitarized BFKL Pomeron as solution, in the LL(l/x) approximation. 
The scattering of the onium qq (dipole) with the nucleus in the rest frame, takes place 
through a cascade of soft gluons, which once taken in the Nc ^ oo limit is simplified by 
the suppression of non-planar diagrams. The gluons are replaced by qq pairs and the dipole 
Mueller's technique for the perturbative cascade can be employed ||3^ . 



The Kovchegov formulation, as the AGL, is a perturbative QCD calculation and the 
considered dipoles from the cascade interact independently with the nucleus. The onium- 
onium frontal scattering has the cross section a = ~2XmA, where the amplitude 

A= —i d'^x dz d^Xi / dzi $(a;, z) F $(a^, zi) , 

where $(x, z) is the square of the onium wave function, x is the transverse separation of the 
qq pair, and z is the longitudinal fraction of momentum of the quark. For the exchange of 
only two gluons, without gluon ladder evolution the function F is [^ 

F(°)(f, xl) = - ""'^^;~^^ x^< (1 + ln(:^)) , (30) 

where x>(x<) is the biggest (smaller) between |x| and \xi\. The two gluons approximation 
is energy independent, but for high energy the contributions of order (a^y)" should be 
included {Y = Ins/M^ is the rapidity and M is the onium mass), since they generate the 
perturbative cascade evolution. The dipole approximation introduces an arbitrary number 
of soft gluons in the square of the onia wave function $, and keeping F as an exchange of 
two gluons avoids to deal with the reggeization of the gluons and the effective vertex. The 
transverse coordinates of the quark and antiquark of an ultrarelativistic onium state in + 
direction are xq = and x, and successively in the evolution the next emitted gluon should 
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be softer. We have p — k and k, as the momenta for the pair, and zi = ^ (in hght-cone 



variables [PH]), having as wave function 



^(°)(xoi, 



^1 



£ki 

(2^ 



gifc;xoi^(o)(^^^^^> 



(31) 



where Xqi = xl — Xq, and ^'■''^ = |\l/'^'')p, keeping factorization in this procedure. This allows 
to obtain the dipole density, n, considering xo2 > P, 3^12 > P, where p is an ultraviolet cut 
also implied by C in the expression below 



n (xoi, x, Y) = x6{x — Xqi) exp 



+ 



a.Nr, 



01 X'2 



2 / 2 2/ 

TT JC Xo2 3;]^2 "'0 



Y 



dyexp 



_2a^^ x^^ 

TT p 

[Y -y) In , 

TT p 



X n(xi2,x,y) 



for y = Ins/M^, which is represented in Fig. (^. 




(32) 



FIG. 7. 



The next step is to obtain an evolution equation for the dipole density assuming the 
propagation of the dipoles in the target is represented by the function 71(0;, 6), where 
h is the impact parameter, and which should be added to the density n2, equally con- 
voluted with 72(3^, &i, 3^2, ^2)5 etc. Assuming no correlation among the dipoles 7n(---) = 
7i(xi, 6i)...7i(a^, 6„), the cross section for the interaction onium nucleus N^x'^i^h^.Y) is 
then given by [^ 



- N{xm,hQ,Y) = XI / ni{xQi,Y,hi,xi,...,h 



, Xi 



X 



2'KXi 



lixi, hi 



\Aj tXj 1} 



■>■ J2 



2-nxj 



d%. 



(33) 



Finally, omitting some steps of the calculation W^ the evolution equation for 
N{xoi,bo,Y) is 



23 



N{xoi,bo, Y) = -7(a;Ki, &o) exp 



4asCF ^,Xoi; 



TT 



P 



+ 



a,,C 



^2 



Y 



dyexp 



ln( — ){Y -y) 

TT p 

2N{xo2,bo + -xi2,y) - 
- N{xo2, bo + -xu, y)N{xu, bo - -X20, y) 



X / £X2 "'"^ 



•^02-^12 



(34) 



where Xij = Xi — Xj, the size of the dipole whose quark has transverse coordinate Xi, and 
the antiquark Xj, 7(0:01, &o) is the propagator of the pair qq through the nucleus, describing 
the muhiple rescattering of the dipole with the nucleons within the nucleus. We denote this 
equation as the K equation. 

The physical representation is comparable with the approach Glauber-Mueller since the 
incident photon generates a qq that subsequently emits a gluon cascade further interacting 
with the nucleus. At large Nc limit the gluon can be represented by a qq pair, and we 
can expect in this limit and DLA that the gluon cascade could be interpreted as a dipole 
cascade. Although beguinning the formulations with distinct degrees of freedom both K 
and AGL resum the multiple rescatterings in their respectives degrees of freedom, which 
allows to consider they should coincide in a suitable common kinematical limit, which we 
will show later on. 

In DLA, where the photon scale of momentum Q^ is bigger than Aq^^, the K equation 
simplifies to 



dN{xoi,bo,Y) _ asCp 2 



dY 



Tl 



01 



"^/^QCD d?X02 



\x, 



2 \2 
02 J 



X 



2 N{x-o2, bo, Y) - N{xo2, bo, Y) iV(xK2, &o, Y) 



(35) 



which is the evolution in transverse size of the dipoles from xqi up to 1/Aqcd- Now deriving 
in ln(l/xoi Aq^^) results 

d''N{x-oubo,Y) asCp 



dYd\n{l/xl,AlcD) ^ 



2-N{xoi,bo,Y) 
xN{xoi,bo,Y). 



(36) 
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setting that the successive emission of dipoles generates larger transverse size for each higher 
generation. 

The hnear term reproduces BFKL at low density, and the quadratic term introduces 
UC unitaryzing the BFKL Pomeron and the equation reproduces GLR once we assume N 
directly related to the gluon distribution function. 

D. The MV-JKWL Formulation 

In the MV-JKWL formulation [0,^ a very dense system is treated in the light-cone 
and considering the light-cone gauge (^4+ = 0), a; = iGiuonl Q%udeon- I^ -^^f- [0 ^^^ gluons 
distribution for small x is proposed for a large nucleus where the degrees of freedom are 
virtual quanta from a classical field generated by the color charge of the valence partons 
(static sources). The approach is originally non-perturbative and the nucleus is considered 
in the infinite momentum frame, transfering the scale of the problem to A = l/vri?^ dN/dy, 
where N is the density of gluons. For small x and a large nucleus ^^(A) is small allowing 
some perturbative calculation in this effective lagrangian formulation for gluons condensates. 

The density of gluons in momentum space is obtained in terms of the correlation of the 
gluons fields, in the light cone gauge. The intrinsic quantum fluctuations are replaced by a 
classical average on the color charge ensemble. The gluons distributions at a given virtuality 
Q^ and x is obtained from the density of gluons in the momenta space dN/dq^q'^q, which is 
a function of the gluon condensate < A°;{x~ ,x) A{x', x') >, being 

xg{x,Q'^)=J d^QX-i — ^- (37) 

The gluons distributions, in this framework where a large number of color charges gener- 
ates a QCD vector potential, is obtained in lowest order by solving the classical Yang-Mills 
equations, D^F^" = j". 

Introducing a regulator in the valence partons current singularity by considering the 
color density p a function of rapidity, it was obtained PUI an analytical solution for the 
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classical correlations, with the property that for high transverse momentum the classical 
gluons distribution obeys the Weizsacker- Williams form, and has its behavior softened as 
lia{k^ / x{y , k"^)) . Here x = !^ l^'^iv^Q'^) is the squared color charge per unity of area for 
rapidity bigger than y. 

In the classical MV the non-linear effects are included in the charge density p solution 
of Y-M equations. The quantum corrections are to be considered, and from Ref. ||3^ the 
perturbative result for the gluons distribution up to second order in as is given by 



1 dN 
•nB? dxd?kf 



a.,r 



1 



xkl 






(38) 



where V = ^ l~ and fx^ is the square of the color charge average density (per unity of 
area). The additional effect of including the hard gluon was treated in Ref. p^, resulting 
in the low density limit the BFKL equation, and for high virtualities the DGLAP equation. 



For high density a complete solution was not yet obtained. In Ref. [|^ JKLW analysed their 
evolution equation in DLA obtaining a generalization of GLR. 

As a summary of the formulations for hdQCD at present, the Fig. (§) presents their 
different regions of applicability as fas as k, in concerned in the ln(l/x) versus Q^ plane. 
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The main questions at this point can be: 

• which is the most suitable form to introduce the UC ? 

• can we relate the distinct formulations in a common limit analytically ? 

• what do we look at the observables as a signature for the UC ? 

The last two questions, I will briefly address in the rest of this presentation following our 
personal contributions to this investigation. 

V. CONNECTION AMONG THE FORMULATIONS 

The AGL equation was originally obtained from the Glauber-Mueller approach, but it 
can be also derived from the dipole representation |^T|. We obtained the cross section for 
the virtual probe G* with the nucleus o"*^*"^ = /q dz j ^-ii|^^*pcr99+^^ that can be expressed 
by means of the dipoles qq once we remind a^^+A _ (^(7^/(7^)cr'?9+^. Now in order to 
estimate the UC the rescatterings of the qq pair into the nucleus should be considered, 
having in mind that |^ cr^ = ^(3as(4/r^)/4)7r^r^X(?Ar(x, 4/r^), where xqn is the nucleon 



gluon distribution. The wave function \E''^* calculated in p6| , ^ is such that I^E'^'P = 
^^^[{e^K^^ett) - eKi{ert)/rtf + l/rt{eK^{ert))% where e" = Q^z{l - z), and the K, are 
the modified Bessel functions. For small z and er^ << 1 we obtain 

4 Ca f^ dx' f°^ d'^n 



xgA{x,Q ) = , , / , ^ 4 



X 2 



1 _ e-i-?.'^(^*) . (39) 



From this equation we can obtain the AGL equation in the dipole representation by 
differentiating in y = In 1/x and InQ'^/Aqf-,^ having 



dydlnQ^KlcD J ^ 

valid in DLA, considering each gluon of the cascade as a qq in the high N^ limit, and where 
C = 2Ca/tt^Cf. For a central collision and S± = ttR^ and 5(0) = A/ttR^ 
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d'^xg{x, Q"^ 



DQ' 



\ — e JVcSj^Q 



-1^9 A 



(41) 



dydlnQyAlcD 
for Nc = 3,Cf = NJ2 at high A^^, and where D = ^^^f^' . 

The GLR for a cyhndrical nucleus is immediately obtained from Eq. (|41| ) by its expansion 
up to second order in xqa- Again, for small UC only the first term contributes which 
reproduces DGLAP in the DLA limit. Those results are in Ref. El 



Now, the Eq. (34) is the K equation that in DLA, where the scale of momentum of the 
photon Q^ is higher than the scale of momentum of the nucleus Aq^^,, simplifies as 






dY 



TC 



^01 



dXg2 



,2 
01 



[X, 



2 \2 
02 J 



X 



2N{x-o2,bo,Y) - N\x-o2,bo,Y) 



(42) 



This equation considers the evolution of the dipoles from xqi up to 1/Aqcd in the 
transverse direction. Now deriving Eq. (^) in \ia{l/{xl-^^AQ^^)) we get 

d^Nixoi,bo,Y) a.Cf 



Xs'^F 



N]N. 



(43) 



dYdHl/xl.Alcn) ^ ' 

We should relate now the function N{xoi,bQ,Y) with the gluon distribution function. 
For that we consider the structure function F2 for the nucleus as obtained in [O], following 
25 , and analized in [^ for bt = 0, which is 



F2^{x,Q' 



^" R^fdzf^m^ 



An'^a^ 



IT 



x2 






(44) 



This estimates the UC for the nuclear structure function, for central collisions in the 
DLA limit in the Glauber-Mueller approach. 

Considering the unitarity corrections due to the multiple rescattering of the qq pairs with 
the distinct nucleons into the nucleus, from the just obtained expression for F^, it results 
the relation 



Nixoi,bo = 0,Y) 



O c _L 



(45) 
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where xqi = xq — xi = rt, Y = ln(s/Q^) = ln(l/x) establishing a connection between the 
cross-section of the qq pair and the gluon structure function in DLA hmit. 

We are in good terms to verify the connection among the K and AGL formulations since 
we already obtained the AGL equation in the dipole formulation, Eq. (40), the K equation 
in the DLA limit, Eq. (42), the cross section of the pair through the dipole density from K 
and the nuclear gluon distribution function, Eq. (45). 

Having Eq. (45) in Eq. (41) and for xqi ~ '2/Q, as in [l^, we obtain 

d'^xgA{x,Q^) 



DQ' 



I — e ^csxQ 



7^9" 



(46) 



dydHQyAlcD) 
result already obtained, and that gives GLR as a limit. 

Our comparison has physical meaning for dipoles with small transverse sizes and for the 
above connection among A^ and xqa [Eq. (45)]. 

In Refs. |]38| , |39[| it was applied the Wilson renormalization group to the model of McLerran 
and Venugopalan. The non-linear evolution equation then obtained deals with the weight 
function of the color charge densities, valid at leading order ag and for densities up to l/a^. 
The complete analytical solution is not yet obtained but some limits are discussed. At low 
densities BFKL is recovered, and then at DLA at large Q^ DGLAP is recovered. In the 
work [^ is proposed the equation 



d'^xg{x,Q'^,bt 
dyde 



PQ' 



1 



X 



exp{l/K)Ei{l/K) 



(47) 



where /? = N^{N^ - l)/2 and k{x, Q^, h) = 2ajTi{N^ - l)Q'^xg{x, Q^ h). 

For large k, a factorized ht dependence and considering a central collision we obtain for 
this equation 

d'^xg{x,Q'^^ 



dyde 



PR'Q\ 



(48) 



which solution is 



xg{x,Q^) = (37rR^Q^\n{l/x) , 



(49) 
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presenting the same Q^ and x behavior as the asymptotic solution for AGL. The main point 
is the partial saturation of the gluon distribution presented in both formulations in the 
asymptotic region. A connection among those two formulations in a more broad kinematical 
region is still an open question. 

The asymptotic behavior of the structure function also required our attention. Consid- 



ering the relation of cr'''^ and xg[x^ Q ) we can write ||29 



F2{X.Q') 



2a, fQ' dQ^^ 



xg{x,Q'^ 



(50) 



9tc JQi Q2 

which is a leading twist equation, with limited application for high densities, due to higher 
twist terms related with the UC 

Using the solution of AGL in the asymptotic regime as input in the above equation 
we obtain F2{x,Q'^) ~ ^R'^Q'^ln{l/x), which again presents partial saturation, meaning 



the Froissart limit is not violated |Q. Analogous result was obtained by Kovchegov |^ 



employing the solution of the K equation ||TI[ . We obtained that the asymptotic behavior 
of F2 is a general characteristic that appears to be independent from the approach that is 
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Assuming the asymptotic behavior of the gluon function is xg{x, Q^) = 2Q'^R'^/3TTas, it 
implies saturation for F2 (~ R^Q^) for very small x. However this result should be taken 
with caution since it is valid in a kinematical region where higher order in the partonic 
density are not significative. The subject of saturation is a tricky one and it seems we are 
far from establishing its features in a solid theoretical basis [Q. Important contributions 
to these challenging aspects of hdQGD are to be found in Mueller EM for the theoretical 



discussion and Golec-Biernat and Wiisthoff |24| for a phenomeno logical application. 



In ||43|j44l we were able to show that 



F^M') = ^.i: 



e?/''^°^[C + ln., + E,(., 



1/Q2 Tir^ 



qi\ 1 



(51) 



where Kq = A/dHg. From that we can estimate the UG for F2 in the DLA limit. For large 



Kg, and using the asymptotic solution of AGL, we obtained [^ 
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37r2 



F2(x,Q2)^^ln[:^ln(l/x)], (52) 



when higher twist terms are considered in F2. This is a softer behavior, but in both cases 
there is no violation of the Froissart hmit. This above equation was not studied in K or 
MV-JKLW approaches. 

From the aheady obtained results it follows the identity [^ 

d\nQ^ -^2(x,g), (53) 

as an important signature of the asymptotic regime of QCD for dense systems. It is relevant 
to mention that for the same center of mass energy this regime is reached for nucleus for 
smaller partonic densities than in the nucleon case, since ka = A^/'^k^. 

VI. PHENOMENOLOGY 

From the Glauber-Mueller formalism for high dense partonic systems was demonstrated 
the AGL equation and its asymptotic behavior. It was also obtained the nucleon and the 
nuclear gluon distribution function as well as the respective structure functions and deriva- 
tives. This formulation incorporates the UC required by the Froissart bound, through a non 
linear dynamics. 

In this section the behavior of the main observables obtained in ep collisions, and relevant 
for eA collisions, will be analysed with the goal to shed some light in the subject of UC. 

For ep we studied the behavior of the proton structure function F2, its derivative -r^^hn 
the charmed component of the structure function F|, and the longitudinal distribution 
function, F^ P5|. 



There is a large amount of data from HERA to motivate a detailed study of these observ- 
ables directly connected with the gluon distribution function. As previously demonstrated 
the gluon distribution is modified in a unitarity corrected formulation, meaning those ob- 
servables should be affected. 
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We were lucky to show that the qi^q2 , the F2 are clearly modified. Also, the eA analysis 
provides stricking results for the nuclear structure function F^ and its derivative, as an 
important signature of the UC corrections. These results are important since they are a 
prediction both for HERA-A and for e-RHIC, in which a high dense parton system should 
be formed. 

The increasing of F2 in HERA in the small x region (10^^ > x > 10^^) is observed even 
for small virtualities {Q"^ ~ 1 GeV^). Taking F2 ~ x~'^, for small x data is compatible with 
A = 0.15 (g2 = 0.85 GeV^) up to A = 0.4 (Q^ = 20 GeV^). This is described by DGLAP 



with suitable input initial condition for Q^ and distributions by different groups [^,^]. It 
will conduct to the idea UC are not observable in the HERA kinematical range. We have 
shown that the structure function is too inclusive in the gluon function to clearly explicitate 
the UC. 

We arrived at a different conclusion applying AGL to offfe, -Pl and F|, all observables 
directly associated with the gluon function. 

The derivative of F2 is 

dF2 R^Q"^ 



aing2 27r2 



5:e2[C + lnK)+EiK)], (54) 



that we solved using the same procedure as Eq. (20). The usual parametrizations p7|j47[ 
do not include the UC for the gluons explicitly. We use Eq. (20) for A = 1 as input, and we 
obtain the corrections from both sectors quark and gluons. The last one gains in importance 
as Q^ increases. In Fig. (||) the results for ^f^ are presented for R^ = h GeV^. For the 



complete discussion we refer to |^8|. The UC for both sectors are able to describe properly 



the data including the turn-over. Our conclusions is this is a good observable to evidentiate 



the presence of the UC. This question was addressed also in ^^ calculating the suppression 
factors separately. 
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FIG. 9. 

We believe that the UC should be extracted from data related to observables that are 
directly dependent of the gluon function. The longitudinal structure function Fl is a difficult 
measurement requiring distinct values of the center of mass energy, meaning different energy 
beams. An alternative is to consider the radiation of a hard photon from the incident 
electron, reducing the center of mass energy. If this should be done we can study -Fl(x, Q"^) 



in the small x region |43 



Expressed considering the quarks transverse momenta due to gluon radiation, the longi- 
tudinal structure function reads 



Fl{x,Q' 



2tt 



-X 



X y^ 

X, 



lF,{y,Q') 






(55) 



where y = Q'^/sx and the dependence on the gluon distribution is explicit, meaning this 
function should be sensitive to unitarity corrections in HERA kinematical region. Our results 
for small x region are in Fig. ([ToD ||43|] compared with the HI data [^. Although it seems 
to be a good observable to evidentiate the UC the available data do not allow any definite 
conclusion for now. 



33 



Fl 1.0 




0.00010 0.00020 0.00030 0.00040 0.00050 



FIG. 10. 



A problably more promissing observable is the rate Rp = -^2^(2;, Q'^)/F2{x, Q"^), where F2 
is the charmed component of the structure function. Considering the approach of boson- 
gluon in order to create the cc pair we obtained in the Glauber-Mueller formalism the ratio 
Rp. 



This ratio is presented in Fig. (11) |43| as a function of ln(l/x). There is strong modifi- 
cation of the ratio once UC are included in the calculation. We urge data in this observable. 
The suppression is much stronger than in the F2 case, and we expect a lower production of 
quark charm for small x, and this is related with the production of J/\l/ which is proportional 
to the square of the gluon function. 
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FIG. 11. 



Finally, one of our most stricking results concerns eA physics, and is related with the 
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high dense partonic system in the nuclear medium. The nuclear shadowing is a challenge for 
hdQCD and mainly important for HERA-A, RHIC and LHC physics. We estimated how 
the nuclear structure function and its derivative are modified by the effects of high partonic 
densitiy. 

The shadowing corrections to F2 are associated to the rescatterings of the qq in the 
nucleons into the nucleus, being dependent on the nucleon gluon distribution function. Here 
also we separate the two cases: quark sector, where the gluon distribution is not modified by 
UC, and quark + gluon sector, where now the gluon distribution is modified a la Glauber- 
Mueller. The results are presented in Fig. (|12D ^^ as a ratio Ri = F^/AF<^ , showing that 
for small x the gluon sector contribution should be included, and promote saturation. 
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FIG. 12. 

We obtain that the suppression due to the shadowing in F^ is proportionally smaller 
than in xqa in a perturbative framework, in a different result that in [|l^ where soft physics 
is the main issue. As a new result the saturation of the ratio is attained at HERA-A region 
when the gluon sector is included. The presence of saturation in the perturbative region 
{Q^ > 1 GeV^) denotes the large shadowing corrections in the gluon sector. 

The analysis was extended to the derivative of the nuclear structure function 






(56) 



considering the contributions of the quark and the gluon sectors to the UC. for HERA-A 
s = 9.W GeV2. 
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The predictions are in Fig. ([T3|) [p^ compared with a DGLAP calculation with GRV 
without nuclear effect. The expected turn-over is present in the orthodox calculation but it 
is A independent. The behavior of the derivative is different once UC are considered since 
the maximum is A dependent and runs to higher values of x and Q^ as A increases. We 
conclude this is the best quantity to look for unitarity corrections, evidentiating the same 
partonic density is reached as A increases for higher values of x and higher values of Q^, 
corroborating a perturbative calculation. 

This is a strong motivation to develop this calculation for heavy ion physics, and try to 
connect this formulation with the research in AA physics where the quark-gluon plasma is 
expected to be produced. 
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VII. OUTLOOK 



Several aspects of the formulations for hdQCD and in our approach to the subject as 
well require further investigation. I understand the formulation of high dense partonic 
system should incorporate the methods of nonperturbative physics and non-linear dynamics 
in order to present a comprehensive formulation for a large kinematical regime in x and Q^, 
besides incorporating the A dependence. However significative progress in the description 
of hdQCD has been made in the recent years towards a unified theoretical framework. 
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Particularlly relevant is the role of initial conditions for UC for the perturbative treatments 
and the determination of saturation region, Q^, still to be obtained analytically. Also a 
complete solution of the generalized evolution equation (Eq. (26)) for ^^(Q^) outside the 
asymptotic region is not available. Reaching these goals will allow us to have a more complete 
dynamical description of the non-linear phenomena of transition between large distance and 
short distance physics promoting QCD to a more understandable and applicable theory. 
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